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Abstract 


We derive the exact Helmholtz free energy (HFE) of the standard and 
staggered one-dimensional Blume-Emery-Griffiths (BEG) model in the presence 
of an external longitudinal magnetic field. We discuss in detail the thermo¬ 
dynamic behavior of the ferromagnetic version of the model, which exhibits 
magnetic field-dependent plateaux in the ^-component of its magnetization at 
low temperatures. We also study the behavior of its specific heat and en¬ 
tropy, both per site, at finite temperature. The degeneracy of the ground 
state, at T = 0, along the lines that separate distinct phases in the phase 
diagram of the ferromagnetic BEG model is calculated, extending the study 
of the phase diagram of the spin-1 antiferromagnetic (AF) Ising model in 
[S.M. de Souza and M.T. Thomaz, J. of Mag. and Mag. Mat. 354 (2014) 205]. 
We explore the implications of the equality of phase diagrams, at T = 0, of the 
ferromagnetic BEG model with ^ = — 2 and of the spin-1 AF Ising model for 
D 1 

PT > 2- 
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1 Introduction 


For a long time simple 1-D spin models have been used as toy models for a better un¬ 
derstanding of real systems with coupled spins. Experimental verification of the results 
derived from those toy spin models is difficult, given the complexity of real spin systems 
for any spatial dimension. The development of optical devices permitted the simulation 
of a few 1-D spin models in arrays of cold atoms. In 2011 Simon et a/.|Ij simulated 
the 1-D spin-1/2 Ising model in the presence of a magnetic held with longitudinal and 
transverse components at low temperature. Such possibility encourage us to explore the 
thermodynamic characteristics of 1-D models. 

Recently one of the authors applied the transfer matrix methodJ2J [3, 0j to the calcula¬ 
tion of the exact thermodynamics of the 1-D spin-1 Ising model, with single-ion anisotropy 
term, in the presence of an external longitudinal magnetic field[5]. The present work ex¬ 
tends that discussion to the classical 1-D Blume-Emery-Griffiths (BEG) model [6] with 
external longitudinal magnetic held. This model is classical and its exact thermodynam¬ 
ics can also be derived by the transfer matrix method. The presence of an extra term with 
respect to the 1-D Ising model with single-ion anisotropy term can modify the behavior 
of the quantum chain, mainly its phase diagram at T = 0. In the present article we study 
the thermodynamics of the one-dimensional BEG model in the presence of an external 
longitudinal magnetic held. The phase diagram of the model is discussed in detail for the 
ferromagnetic case, and for two different regions of the parameter pj, and complemented 
by the discussion on the phase diagram of the spin-1 AF Ising model[5], 

In section [2] we present the Hamiltonian of the standard BEG model in the presence of 
an external longitudinal magnetic held. We show the relation between the Hamiltonians 
of the standard and staggered versions of this model, to be used in relating their thermo¬ 
dynamics. In section [3] we discuss the phase diagram, at T = 0, of the ferromagnetic BEG 
model. Its thermodynamics is presented in subsection 13 .1 1 through the behavior of three 
thermodynamic functions per site: the ^-component of the magnetization, the specihc 
heat and the entropy. The entropy per site along each line that separates distinct phases 
in the diagram of the ferromagnetic BEG model is calculated at T = 0. In section [T] we 
compare the three previous thermodynamic functions of the ferromagnetic BEG model 
with p| = —2 and the spin-1 AF Ising model at very low temperature. We also extend 
the discussion on the phase diagram of the spin-1 AF model in Ref. [5j in order to include 
the degeneracy of the ground state of the model at T — 0. Our conclusions are presented 
in section 0 In Appendix [A] we present the main steps to calculate the exact Helmholtz 
free energy (HFE) of the standard and staggered one-dimensional BEG model for arbi¬ 
trary values of the parameters. The states and energies of the dimers present in the spin 
configurations of the chain are shown in Appendix [B] In Appendix Owe have the ground 
states of the BEG model in the presence of an external longitudinal magnetic held and 
their respective energies. Finally in Appendix [D] we show how to calculate the degeneracy 
of the ground state along the lines that separate the different phases along the diagrams 
of the chain models at T = 0. 
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2 The Hamiltonian and HFE of the 1-D Blume- 
Emery-Griffiths model with a longitudinal mag¬ 
netic field 

Eq.(5) of Ref. [5], 

N 

H s r\J,h,D)=Y J [JS‘St +1 -hSt-hSU 1 + D(St) 2 + D(S; +1 f], (1) 

i =1 

is the Hamiltonian of the one-dimensional classical spin-1 Ising model with the single-ion 
anisotropy term with the crystal held D, the Blnme[7j -CapeljSj model, in the presence of 
an external longitudinal (z-axis) magnetic held h, symmetrized in the nearest neighbours. 
Here, Sf is the ^-component of the spin-1 operator in the i-th site (\S\ 2 = 2), and J is the 
exchange strength. For J > 0 we have the anti-ferromagnetic (AF) version of the model, 
whereas for J < 0 ferromagnetic version is obtained. We assume that h > 0, that D may 
have any real value, and that the chain has N sites and it is periodic, i.e. S^ +l = S*. In 
this paper we use natural units e = m = h = 1. 

Adding the term — iL(S'?) 2 (S' l z +1 ) 2 to Hamiltonian (JTD , with K £ M, yields the Hamil¬ 
tonian of the Blume-Emery-Grifhths (BEG) model[6, [9j 

N 

H BEG (J,h,D,K)=J2 [JStSUi - hS*-hS* +1 + D(S’) 2 + D(S’ +1 ) 2 

i=l 

- K(S?) 2 (S? + J 2 ]. (2) 

This Hamiltonian also satishes the periodic condition. (The Hamiltonian (1) with L — 0 
in Ref. [5] describes the BEG model in a non-symmetrized form.) 

Ref. [5] discusses at length which quantum state(s) is (are) favored by each term on the 
r.h.s. of Hamiltonian ©, regarding the minimization of energy; we will not repeat this 
discussion here. Let s ? be the eigenvalue of the operator S'?. In ([2]), the term in K will, for 
K > 0, favor the dimer states (i.e., relative to two neighboring spins) in which s\ = ±1, 
independently of their relative orientation (they may be either parallel or anti-parallel). 
On the other hand, for K < 0, the favored dimer states will be those with at least one 
null eigenvalue, i.e., sf = 0. Section [3] will describe how the term in K changes the T = 0 
phase diagrams of the classical ferromagnetic spin-1 Ising models presented in Ref. [5], 
The staggered BEG model in its symmetrized version reads 

N 

H TEc(J.,h.,D„K.) = £ [j.s;s^ 1 -(-i) , h.s;-(-i)» 1 h.s^ 1 

i=l 

+ D,(S‘Y + D,(S* +i y - KYSlfiSl,) 2 } ■ ( 3 ) 

This Hamiltonian will also be subject to the spatial periodic condition. We assume that 
the chain has an even number of sites, so N = 2M, in which M £ N. 

The mapping S'? — > (—1 ) l Sf in Hamiltonian ((2D yields the relation between the stan¬ 
dard and staggered BEG models 

H%% G (Js,h a ,D„K a ) = H BEG (-J,h,D,K)- (4) 
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hence they have the same thermodynamics if J s = — J, h s = h, D s = D and K s = K. 
The ferromagnetic staggered BEG model (J s < 0) has the same thermodynamics as the 
AF standard BEG model (J >0). The AF staggered BEG model (J s > 0) has the same 
behavior as the ferromagnetic standard BEG model (J < 0) at any finite temperature. 

From now on we will restrict our discussion to the thermodynamics of the standard 
Hamiltonian (J2J) of the BEG model. The thermodynamic behavior of the staggered BEG 
models at finite temperature can be obtained from the corresponding standard models by 
using (@D- 

In appendix [A] we show the calculation of the exact expression of the Helmholtz free 
energy (HFE) of the ferromagnetic and AF BEG models in the presence of a longitudinal 
magnetic field by the transfer matrix method[21 [3J, [3j, valid at any finite temperature 
T > 0. In 1975 Krinsky and Furman[TU] calculated this thermodynamic function for 
those BEG models. Our expression of the HFE for non-null external longitudinal magnetic 
field h 0 written differently from that of Ref. [TO]; ours has contributions only from real 

functions of the parameters of Hamiltonian (J2D and of {3 = —— , in which k is Boltzmann’s 

KjI 

constant and T is the absolute temperature in kelvin. Although the results derived in that 
appendix are valid for both the ferromagnetic and the AF BEG models, in the following 
sections of this paper the discussion is restricted to the ferromagnetic case. 


3 The phase diagram of the ferromagnetic BEG 
model at T = 0 

The Hamiltonian ([2]) can be written as the sum of Hamiltonians of dimers on neighboring 
sites i and i + 1, i G {1, 2, • • • , N}. For the dimer composed of the (i, i + 1) sites, we have 

BlSlMh.D.K) = js;s^ - hS‘ - hS‘ +1 + D(S‘) 2 + D(S‘ i+l f 

~ K(S‘f(Sl +l f. ( 5 ) 


The ferromagnetic case corresponds to J < 0. 

Let and s\ be the eigenstate and eigenvalue, respectively, of the z-component of 
the spin operator at i-th site, Sf, so that S*\sf)i = with sf G {—1,0,1}. The 

energy £^ i+ 1 of the dimer (i,i + 1), described by the state \D) iii+ i = |s})i ® |sf +1 )j + i is, 
in units of | J\, 


£■ 2 , 2+1 

TT 




S i S i +1 


\J\ 

h_, 

V\ 


S 'i + S 


D 


m) + j7| ( s f) + ( 


^ + i) 2 


K 

W\ 


(s;)\si + 1 ) 


( 6 ) 


with sf, sf +1 G {0, ±1}, and i G {1, 2, • • • ,N}. All the parameters of the Hamiltonian (J2]) 
are scaled in units of | J\: , yjj and —; correspondingly, the inverse of the temperature 

scales as \J\f3. 
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In appendix [B] we present the nine possible dimer configurations of neighbouring sites 
in the chain and their respective energy per unit of \J\. The ground state of the ferro¬ 
magnetic BEG model is composed of dimer states which minimize the energy at T = 0. 

K 

The value of the parameter —— determines the general structure of the T = 0 phase 

Ml 

diagram of the ferromagnetic BEG model. 

K 

i ) The case — < —1. 

Ml 

The T = 0 phase diagram for this case is represented in FigjT^,. It resembles the 
phase diagram of the classical spin-1 AF Ising model with single ion anisotropy term and 
external longitudinal magnetic field, discussed in Ref.[5j. However, the AF Neel states 
|T 0 )g and | v h 0 )G' [cf. Eqs. (IC.3a|) and (IC.3b|) of [J5]] are naturally absent from FigJTM 
The ray (half-line) 171 = 171 extending from the origin separates the phases A and E/E', 


_ D_ 

HI HI 

h 

HI 


whereas the parallel ray yj\ — yyf — 1 — - 777 , extending from the point ( f^-, 0 ), in which 


K_ 

HI 


-rjr = 1 T 777 < 0, separates the phases E/E' and B. Correspondingly, the ray = — A) 


D X 


HI : 


HI 


HI 


HI 


HI 

K 


from the origin separates the phases A and F/F' and the parallel ray ^fi = — TTf + l + fTi 


HI 


HI 


from ^ jtj, 0J separates the phases F/F' and C. The two rays from the origin and their 
parallel rays are displaced horizontally by jjp Such displacement increases as the value 
of j 7 | decreases. 

The chain states corresponding to each phase are described in the Appendix O of 
this paper. Those states are represented as |T 0 )a 7 in which X e {A, B,C, E, E', F, F'}. 
The phases {A,B,C} correspond to the nondegenerate chain states |\&o)aj | v I / o)b and 
|T 0 )g, respectively [cf. Eqs. (IG.laj) (IG. lc[) ] . whereas the phases {E/E ', F/F'} are twofold 
degenerate, corresponding to the pairs of chain states (|T 0 )£;, (To)^) and (|To)f, |T 0 )f')> 
respectively [cf. Eqs. (lG.ld|) -( ]C.lgD . 

The number of degenerate ground states corresponding to the lines and critical points 
separating the phases in FigOJi at T = 0 can also be calculated. Along one such line, we 
may determine the possible states of neighboring sites in the chain, which are those dimer 
states listed in Appendix |B] that minimize the energy along that line. The same guideline 
can be applied to the critical points 7 Z and T■ In appendix [D] the degeneracy of states on 
those lines and points for T = 0 is detailed. 

Let Qu^v(T = 0) be the number of degenerate ground state vectors along the line 
that separates two distinct phases U and V of the phase diagram in FigdM excluding the 
critical points 1Z and T. These degenerate states do not necessarily satisfy the periodic 
spatial boundary condition. Our aim is calculating the entropy per site, at T = 0, along 
the lines that separate the distinct phases in the phase diagrams of FigsJIk and [T]o. In 
Appendix [D] we show that the entropy per site along those lines is the same whether we 
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take into account the periodic spatial boundary condition or not. We obtain 


n B ^c(T = o) 

Qe/e't±f/f'(T = 0 ) 
Qat±e/e'(T = 0 ) 


= 2, (7a) 

= | 2^ (1 + (-1)"+') + 2? (1 + (-1)") . (7b) 

= £Ia^±f/f’(T = 0 ) = Qb^±e/e'(T = 0 ) = VLc^±f/f’(T = 0 ) 


1 

10 


’l + y/S' 


N 


(5 + 3V5) + - 


1 - 


N 


(5-3^5), 


(7c) 


in which N is the total number of sites in the chain. We are assuming that N is even: 
N = 2 M, in which Me N. Moreover, the number of degenerate ground state vectors at 
T — 0, at the critical points 7 Z and T is given by 


Qn(T = 0) = 3 x 2 n 1 , and (8a) 

Q t (T = 0) = I (2 JV+2 - (-1)^) - (8b) 


Again fl^ and fir are the total number of degenerate ground states at these critical points, 
including the states that do not satisfy the periodic spatial boundary condition. 

Except for the line separating the phases B and C in the phase diagram of FigOJi, 
all other lines are highly degenerate. The results (ITaj) - (l8b]) do not depend on the par¬ 
ticular value and sign of the exchange strength J; rather, they depend only on which 
dimer configurations yield the minimum energy for the parameter scenario along each 
line. Consequently, a comparison of the phase diagram in FigOJi and the phase diagram 
shown in Fig. lb of Ref. [5] for the AF spin-1 model with single ion-anisotropy and external 
longitudinal magnetic held shows that the phases and the lines separating them are the 
same. Moreover, those lines have the same degeneracy in both diagrams. 


K 

ii) The case — > — 1. 

|d| 

The T = 0 phase diagram of the ferromagnetic BEG model in this case is shown in 
FigHJo. The lines jj-\ = ^ J 7 f ~F 2 | 7 [, separate the phases A and B, and the phases A 

and C , respectively. All phases A, B and C are nondegenerate, and they are described 
by the chain state vectors |^ 0 )^, |T 0 ) B and |T 0 )c, respectively. 

Appendix [B] shows the nine possible dimer states and their corresponding energies. 
Along the horizontal line that separates the phases B and C in Figjl}) to the point 

the energies of the dimers \D^) iji+1 and 1-0^)^+! [cf. (IB. 2 II and (IB.311 ] are the same 
and correspond to the minimum value of energy out of the nine possibilities. The chain 
ground states that can be constructed from \D^)i ji+ i and \D { - C ' ) ) i)i+ i along that line are 
|To)_b and |To)c [cf- (IC.lbjl and (IC.lcjl of Appendix IC], The degree of degeneracy of the 
ground state along this line is equal to 2. 

Likewise, along the line separating the phases A and B in Fig|T}y the dimer states 
with minimum energy are |72^)j ii+1 and |[cf. (IB.II) and (IB.21) ]. and from those 
the chain ground states correspond to |^7 0 )a and |To)_b [cf. (IC.laj) and (IC.lbj) . Hence the 
degree of degeneracy of the ground state along this line is also equal to 2. 
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A similar argument holds for the line separating the phases A and C ; its degeneracy 
is also equal to 2. 

The existence or not of an exponentially growing degeneracy of the ground states along 
the separation lines in the phase diagrams of FigsJIb and db, at T — 0, determines the 
thermodynamic behavior of the ferromagnetic BEG model. 


3.1 Thermodynamic behavior of the ferromagnetic BEG model 

In this subsection we discuss three thermodynamic functions per site of the ferromagnetic 
(J < 0) BEG model: the ^-component of the magnetization, 

1 <9W 

M z {J,h,D,K ] 0) = --- w , (9a) 

the specific heat 

C(J,h,D,K-p) = -p 2 ^p-, (9b) 


and the entropy 


S(J, h, D,K] (3) 


0 


2 dw 

d/3 ’ 


(9c) 


in which W( J, h, D , K] (3) is the HFE of the model. Its exact expression at finite temper¬ 
ature T (in kelvin) can be found in appendix [AJ 

In the following we let J = — 1; the remaining parameters in the Hamiltonian ([2]) are 
in units of |J|: A, yjj and The inverse of temperature (3 is scaled as \J\/3 in those 
functions. Our aim is verifying the effect of the degeneracy along the separating lines in 
the T = 0 phase diagrams of FigsJIb anddb on these thermodynamic functions. We will 
discuss separately the regions of the phase diagram in which m < -1 and m ^ - L 


h_ 

W\ 


for J = — 1 and ^ = 

—1.5. That corresponds to examining, in the phase diagram of FigHb, a vertical straight 
line in the upper half-plane yjy > 0. Two scenarios can be qualitatively distinguished. 


The case ^ < —1 is exemplified in FigJ2I which shows AT 


The first scenario is exemplified by Fig|2f, which shows A4 Z (A 


for 33- = 
° HI 


-2 and 


\J\(3 G {1.5,5,100}. This figure corresponds to following upwards the vertical line Li in 
the phase diagram in Fig{Tb, which begins at the point corresponding to y- — _2 and 
h = 0, on the horizontal axis. Along the vertical line L i, the degree of degeneracy of 
the ground state vector is 2. Figj2b shows that for (J = —1, A = —1.5 and = —2), 

the step-function shape of M. z (] 7 j) is lost about \J\/3 = 1.5. The second scenario is 

exemplified by Fig|2b, which shows M z (jjy) for T = 0.25 and \ J\(3 G {10,100,1000}. 
It corresponds to following upwards the vertical line L 2 in the phase diagram in Fig Jib, 
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starting from the point corresponding to py — 0.25 and h — 0, on the horizontal axis. L 2 
intersects two phase transition lines: from A to E / E' at py = 0.25 and from E/E' to B at 
py = 0.75. The degeneracies of the ground state vectors along these two transition lines 
are exponential [cf. Eq. fjTc]) ]. Consequently, the shape of the curve .M z (py) at \J\/3 = 10 
is quite distinct from that of the corresponding curve at | J\/3 = 1000. For the latter, the 
^-component of the magnetization has two plateaux that resemble quite closely a sequence 
of step-functions. At \J\/3 = 10, however, the plateaux are no longer present. That same 

behavior of A4 Z ^py^ f° r different temperatures is seen in any curve for py > pj-, in which 

pt = 1 + R’ and R <-L 

Figsj3]show the function Af 2 (py) for py > — 1. In Figs]3}i. and |3h> we have J = — 1 
(ferro) and py = 1. In Figj3k we choose py = —0.5. The vertical line in the phase 
diagram of Fig{T}o followed by the argument py of the function M. z , begins at py = —0.5 
with h — 0 and crosses the phase B of this diagram at T = 0. From FigJ2K we see that 
the function _M z (py) has the plateau Ai z = 1 up to \ J\f3 ~ 5. For | J\/3 = 1.5 the curve 
Af 2 (py) in FigJSi differs from the step-function around py ~ 0. One is reminded that the 
degeneracy of the ground state vectors along the line that separates the phases B and C 
in Fig {lb is equal to 2. 

In FigJSJo we have py = 1.8. The function _M 2 (py) is plotted for the variable py varying 
along the vertical line in the phase diagram^ with py = 1.8. This vertical line crosses the 
phase A and cut the line that separates the phases A and B in the phase diagram [TJo at 
py = 0.8, going along the phase B. The number of degenerate ground state vectors along 
the line separating the phases A and B in the phase diagram is equal to 2. In Figj3b the 
curves of the function Af z (py) has two plateaux, A4 Z = 0 and M. z = 1, for \J\/3 up to 5. 

It looses the step-function form for | J\/3 ~ 1.5 (a high temperature T ~ ^). Comparing 
the curves of the ^-component of the magnetization as a function of the py in FigsJ2£i, 
[3]a and EJo, we verify that the plateaux in those curves are present up to a high enough 
temperature. The common fact about these three curves is that the vertical line in the 
phase diagrams of Fig. |Tj(for an increasing value of py and a steady value of py) cut phase 
transition lines for which the degeneracy of the ground state vectors are not exponential. 

The entropy per site 5( J, h, D , K] /3 ) can be derived from the HFE W( J, h, D , K ; f3 ) 
of the model through the relation between these two functions presented previously, see 
eq. (15c]l . or through the number of states with energy between E and E + SE, with SE -C 
F[T3], 12(F), in the thermodynamic limit, 

S(J, h, D, K; f5) = lim ln(n(E)). (10) 

JV-xx) iv 

From this point on, we will use a system of units in which k — 1. 

Using an algebraic manipulation program we derive from the results of appendix [A] 
the temperature dependence of the entropy per site of the BEG model (ferromagnetic 
and AF models), valid for T > 0. In eqs. dTal) - (l8b|) we present the number of degenerate 
ground state vectors along the line that separate the different phases in the diagrams [1^ 
anddb, at T = 0, and the critical points 1Z and 77 


By varying the inverse of temperature, \J\/3, up to 10', we verify numerically that 
the entropy per site, «S(—1, A, jj-r, | J|/3), has a strong indication that the limit of this 
thermodynamic function as T —y 0 (\J\/3 —> oo) along the lines that separate the phases 
in FigsJTJi anddb are as follows. 

1) For yjy < —1 (phase diagram [Qi) we have, for the B C transition line, 


IJeL 5 ^- 1 ' 0 ']7|’|7[ ; I J l« = F" 7f ln ^-o(T = 0)) = 0. 

with ttt < ttt and rh — 1 + - For the E/E' ^ F/F' transition line, 


(11a) 


|J| M 


|J| 


ij|j™« 1 5E/E - F/F '* _1 ’ 0 ’]7f’ j7f ;|J| ^ - jfc, If - o)) 

= i(n(2) ss 0.3466, (lib) 


with -^ < yy| < 0. For the B 


E/E' transition line, 


iJEL S b = e/e ' (_1 ’ |7f 1 ]7f’ pf : W) ~ uli I 3co' Sc=E/E '* _1 ’ J7|’ p|’ pf’ |J| ^ ) 

= lim ln(Q BT±E/E >(T = 0)) 

N —boo iv 

1 


= lim ln(n c ^ F/F ,(T = 0)) 

Af->oo iV 


= In 


'l + 


0.4812, 


(11c) 


with jjj = ±yjf =f 1 q= -pjy and pp > -pp. For the A E/E' transition line, 


-|J| 


I J\ 


\J I ^ kl 


«£.»'pf- pi’ pf; i-pf pi ■ pi : W) 

= lim ln(Vt A ^ E/E i(T = 0)) 

TV—>-oo iv 

= lim i ln(tt AT±F/F i(T = 0)) 

TV—>-oo iv 


= In 


'l +V5 N 


0.4812, 


(lid) 


with A = ±171 and R > °- 

In the multicritical points 1Z and T, at T — 0, we have the limits: 

lim £*(-1,0, yf\,Tyh\J\P) = . f Ji m M" 1 ,0,0, | J\P) 

\J\/3— boo |J| |J| |J|/3—»oo |J| 

= V m Tr = °)) 

N—>oo J\ 

= V m Ji ln[p, T (T = 0)) 

N —boo iv 

= Zn(2) « 0.6931, 


(lie) 
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where p = p = 1 + p (point 1Z) or p = 0 (point T). 

2) For p > 1 (phase diagram [T}r) we also vary the inverse of the temperature, \J\j3, in 
the interval [0,10 s ], and the results strongly indicate that 


lim S(— 1, 

| J\/3 —>-oo 


h D I< 

wvyvi 


; \J\P) = o, 


( 12 ) 


for any value for the parameters p, yy t, when p > — 1. 

Fig J4] represents the entropy per site as a function of jyj- for -pjj = —1.5 at | J\/3 = 10 3 . 

The twin peaks in the curve correspond to the points 1Z (p = —0.5) and T (p = 0) in 
the phase diagram [ljo. Since T > 0 for this picture, the curve is continuous everywhere. 
We verify that at low temperature the points 1Z and T are well featured in the curve 


c f T \ dS(T) 

The relation — = — connec ^ s ^ ie specific heat P er site, i n units of the temper¬ 

ature, and the variation of the entropy per site with the temperature. In what follows we 
will discuss the behavior of the specific heat per site in the two regions of the parameter 

w V\ < -1 and W\ - ~ l - 

For p < — 1 and p < 1 + p, the maximum value of the specific heat per site, as 
function of p at \J\/3 = 10 2 , is of order 10 -85 . For this set of values of the parameters 
p and pr, increasing the variable p corresponds to following a vertical line in the phase 
diagram [T}r that does not cut any line in the diagram that has an exponential degeneracy 
in the ground state at T = 0. The maximum value of C(fTf) \J\P = 10 i s °f order 
10 - ' and at \J\/3 = 2 the highest value of this thermodynamic function is of order 10~ 2 . 
For the temperature varying two orders of magnitude the specific heat per site varies 10 83 
orders of magnitude. 

The specific heat per site with p > — 1 and p e M has the same behavior as described 
in the previous paragraph. The phase diagram of the chain with p > — 1 has no line in 
which the ground state is exponentially degenerate at T = 0. 

In FigIS] the curve C x p is drawn with p = —1.5 and p = 0.25. The variable 
p > 0 follows a vertical line in the phase diagram [1^, at T = 0, that crosses two lines 

that separate the phases A and E / E' at p = 0.25, and E/E' and D at p = 0.75. 
Along these two lines that separate theses phases in diagram [T^l, the ground state vectors 
are exponentially degenerate at T = 0. In FigsJSJr and Eb the curve C(yp) is plotted at 
| J\/3 = 10 3 , but in two intervals of the variable p where this thermodynamic function is 
non null at this temperature. In Fig|5fc the specific heat per site is drawn at | J\/3 = 10. 

We verify from the Figsj5]that the maximum value of the function C(p) at \ J\/3 = 10 3 
is of the same order of magnitude as at \J\f3 = 10. The curve of the specific heat per 
site as a function of p at different temperatures, has the same behavior as described in 
Figsj5]when p < —1 and p > 1 + p. 
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lx AF Ising 


4 Ferromagnetic BEG model at ^ = — 
model at low temperature 

Ref. [5] presents the phase diagram of the spin-1 AF Ising model, with single-ion anisotropy 
term, in the presence of a longitudinal external magnetic field, at T = 0 (see Figj6jr). In 
FigEb we have the phase diagram, also at T = 0, of the ferromagnetic BEG model with 
jj| = —2, in the presence of an external longitudinal magnetic field. Comparing the phase 

diagrams of these two models, we verify that for ^ | the two models have the same 

phases at T = 0. The line l\ in the diagram [6}) corresponds to the vertical line jyj = 
(Having the same phases for yyj > |, however, is not a sufficient condition for the existence 
of a transformation that maps one model onto the other.) 

In subsection 13.11 we discussed the degeneracy of the ground state along the phase 
separating lines in the diagram of the ferromagnetic BEG model, at T — 0, concluding 
that the results (17a]) - d7c|) are independent of the sign of J. A analogous discussion for the 
spin-1 AF Ising model (see FigEb) has not been done in Ref. [5], though. 

The independence of the number of degenerate ground states with respect to the sign 
of the exchange strength J, calculated in subsection 13.11 permits us to affirm that the 
degeneracy of the ground state of the spin-1 AF Ising model along the lines A = with 

jj| > | in diagram Eb, at T = 0, is equal to the result (T7c|) . Again we are also including 
the states that do not satisfy the periodic spatial boundary condition. The degeneracy of 
the ground state of the spin-1 AF Ising model, at T = 0, along the line = ±1 ± 
with yj| > 0, in the phase diagram^ is equal to result (l7b]h under the same situation on 
the spatial boundary condition. 

There are phase transitions at T = 0 in the phase diagram Eb of the spin-1 AF Ising 
model that are absent in the ferromagnetic BEG model with ^ = —2. In the following 
we present the total number of ground states along the lines that separate the phases 
in the spin-1 AF Ising model. These degenerate ground states do not necessarily satisfy 
periodic spatial boundary conditions: 

OZq/g'(T = 0 ) = 

^G^G'^E/E'iT = 0) = 


^b^g/g'(T — 0) — 


The number of degenerate ground states at the critical points (V,V) in the phase 


3, 

Qg^g’^f/f'^T = 0 ) 

N+2 / 1 + ( — 1 )^ 


2“ 

-\AF 


N +1 N-l ( 1 + ( — 1)^ +1 

+ (2—+2—) 1 y ] 


(13a) 
, (13b) 


^C^G/G'(^ — 0) 


1 + 


2-^/5' 




N—l 


+ 1- 
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diagram [ 6)1 is 


n 


AF 

V,V 


2 n +2 


(~l) N+1 

3 


( 14 ) 


We do not show here (the lengthy expression of) the number of degenerate ground 
states at the critical points (Q, Q!) in the phase diagramof the spin-1 AF Ising model, 
at T — 0. Its calculation has been done with the help of a computer algebra system; the 
interested reader may contact the authors for an ASCII hie with that expression. 


The entropy per site of these lines in the phase diagram of the spin-1 AF Ising model 
in Fig j 6 ju , in the thermodynamic limit (.N —y oo), is 


lim 

\J\/3—>oo 


cAF 

°B^G/G 


'( 1 ) 


h D 

W\J\ 


lim 

| J\/3—>oo 


cAF 

°C^G/G 


'( 1 , 


1 K 

= 2 ’Mj 
h D 


— 0; \ J\/3) — 


= In 


l + V^A 

2 J 


« 0.4812 


lim 

| J|/3—>-oo 


Sa^±e/e ’{—1 


h D K 

’PT’PT’PT 


;\m, 


with jj| < ~1 and jyj- > 0. We also have 


(15a) 

(15b) 


Pf-17| • |7[ - 0; I - 

= jjj’ jjj’ jjj = \ J\P) = 2 ^ n (2) ~ 0.3466 (16a) 

= u]!sJ e/e - f,f,{ ~ 1 ’JT\'W\'W W) ' (16b) 

In the phase diagram of the spin-1 AF Ising model we have | < j^j < 1 and 0 < Mr < 

In the phase diagram of the ferromagnetic BEG model, at T — 0, we have ^ < —1. 

The entropy per site of the critical points (V,V), in the thermodynamic limit (N -A 
oo) at T = 0, is 


lim 

| J\/3—>oo 


S'pAF -pi (1, 


h 

V\ 




0; \J\f3) = ln( 2) ^ 0.6931 


lim 

| J\/3 —>-oo 


h D K 

jjj-. |jf> |j| > \ J \P) 


lim S-y( —1, 

|J|/3—>oo 


h 

W\ 



o,4 W ) 


lim 

| J\/3—>oo 


— 1 , 


h 

V\ 


(17a) 


(17b) 
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( 18 ) 


where py = 1 + pr, with p < — 1. We also have at the critical points (Q. Q') 

^ tJj - ±14 = i T?f = 0; | Jl/3) « 0.5886. 

The quantity lim <Sq F q/ can be obtained with arbitrary precision. 

|J|/3—>oo ’ 

The previous information about the degeneracy of the ground state of the spin-1 
Ising model at T = 0 and along the lines that separate the phases in the diagram [6]a, 
complements the information of Ref. j5j about the phase diagram of the spin-1 AF Ising 
model, with single-ion anisotropy term and in the presence of an external longitudinal 
magnetic held, at T = 0. 

Let us examine in what follows the consequences of having the same phase diagram at 
T = 0 for p > | (see the diagrams of FigEJ), in the presence of an external longitudinal 
magnetic held, for both the ferromagnetic BEG model (with p = —2) and the spin-1 AF 
Ising model with single-ion anisotropy term (with p = 0), 

At very low temperature, the contribution to the thermodynamic functions of a model 
comes mainly from its ground state, degenerate or not. In the following we compare 
the ^-component of the magnetization, specihc heat and entropy per site [in eqs. flUall . 
(|9bl) and (EJc]) . respectively] of both models at \J\/3 = 500 and \J\(3 = 100. At the 
latter temperature, we expect that the excited states of the models give contribution to 
the thermodynamic functions due to the exponential degeneracy of these states in each 
model. 

The big difference between the ferromagnetic BEG model, under consideration in this 
section, and the spin-1 AF Ising model is the sign of the exchange strength J. In the 
former (latter) model we have J = —1 (J = 1). We notice from the Hamiltonians ([1]) 
and (J2]) that in order to compare the contribution of the J-term to the partition function 
of the model and the contributions from the external longitudinal magnetic held and the 
single-ion anisotropy term, we have to compare the values of the parameters: J, 2 h and 
2D. 

For a fixed value °f p > |, we have varied the external longitudinal magnetic held in 

whole interval, that is, p G [0, oo) and have compared the three thermodynamic functions 
mentioned previously. 

Along the vertical line in the phase diagrams of FigJHl with hxed value of p , we have 
an interval of p G [0, p ) for which M. z vanishes at \J\j3 = 500 and \J\f3 = 100. In order 
to avoid any singular point in this comparison we dehne the difference of these functions 
as 


- (19) 

FigEJi shows the magnetization difference (fTUD as a function of py for yp = 0.51 and 
\J\/3 = 500, in an interval about A ~ 0.51 with amplitude of order 10~ 10 , in which the 
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transition A ^ E/E' occurs. For p > 0.511, the difference (fT9l) of the functions Ai z is 
a monotonically decreasing function, where we have: 


Diff{M z ) AF (A = 0.52, ^ = 0.51; \J\f5 = 500^ 
Diff(M z ) AF (JL = 1 , ^ = 0.51; \J\P = 500^ 
Diff(M z ) AF = 1.51, = 0.51; \J\f5 = 500^ 


» 4.6 x 10" 14 , 
» 1.6 x 10~ 222 , 
« 3.8 x 10" 444 . 


(20a) 

(20b) 

(20c) 


One is reminded that at p = 1.51 we have the phase transition B ^ E/E' in the 
phase diagrams of FigJHl at T = 0. 

Keeping p = 0.51 and increasing the temperature to \J\/3 = 100, the difference of the 
^-components of the magnetization of the two models decreases; it is shown in FigJ7}:>. For 
Pi ~ 0.51, the value of the parameter p is close to 1 and the ferro- or antiferromagnetic 
nature of the models begins to appear. For p > 0.516, the function Dif f(M z )p F is a 
monotonically decreasing function of the external magnetic field A and 


« 2.5 x 10" 45 , (21a) 

« 9.2 x 10 -90 . (21b) 

For p| = 1 we obtain that the difference (fTUl) of the ^-component of the magnetization 
of the two aforementioned models at \ J\/3 = 500 is \Diff(M~)p F \ < 10 -436 in the interval 
of p G [0.98,1.05] and \Diff(M z ) AF \ < 10 -825 in the interval of p G [1.9,2.1], For 

| J\/3 = 100 these differences decrease and they are < 10 -89 . Increasing the value of p 
(crystal field per unit of |J|) it is impossible to experimentaly measure the difference 
between the z-component of the ferro ( J = —1) and AF (J = 1) models. 

In order to compare the specific heat and the entropy, both per site, of ferromagnetic 
BEG model with p = —2 and the spin-1 AF Ising model, we define the percent differences 

PD(L) "{w W\' w ) = 

LB ™ (-i. pi A A = ~ 2; i'A - Ajk/kl 0; I- 7 '' 3 ) 

LB Fe G (-i,jT\’ww\ = - 2 ’^) 

( 22 ) 



j x 100%, 


Diff(M z ) A / e ( tG = 1,7X7 = 0.51; |J|/3 = 100 


h 


D 


Diff(M z ) AF ( — = 1.51, — = 0.51; | J|/3 = 100 


h 


\J I VI 

D 


\J\ 


\J\ 


where L can be the specihc heat per site, C, or the entropy per site, S. 

First, let us compare the specihc heat per site of the two models. First we consider 
p = 0.51 at \ J\j3 = 500. We obtain that \PD(C)p F \ < 10 -6 % in the interval p G [0, 0.5], 
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but at yyy = 0.51, in which occurs the A ^ E/E' transition in the diagrams of FigjGj at 
T — 0, the value of the specific heat per site at each model is 

C|f G ^-1, A = 0.51, = 0.51, ^ = -2; \J\P = 500^ « 3.04 x 10" 213 (23a) 

C AF 9 (l, |J| = 0.51, = 0.51, A = 0; I J|/3 = 500^) « 8.70 x 10" 8 , (23b) 

showing a difference of 205 orders of magnitude in the specific heat per site of the two 
models. We do not know which mechanism permits so huge difference between these two 
thermodynamic functions. 

For yjy > 0.52, the percent difference (|22|) of the specific heat per site at \J\f3 = 500 is 
\PD{C)j%\ < 10- 8 %. 

Keeping yp = 0.51 and increasing the temperature up to \J\/3 = 100, we plot in Fig|8k 
the specific heat per site versus yp of the ferromagnetic BEG model with yp = —2 and 
the spin-1 AF Ising model in the interval yp G [0.4, 0.6], showing that the two curves do 
not coincide in this whole interval. For yp > 0.8 we have \PD(C)pf\ < 10~ 12 %, and at 
yp = 1.51, when we have the transition of phases B E/E' in diagrams H at T = 0, we 

obtain PD(C )$f = 1.51, -g- = 0.51; | J\fi = 100) ~ 4.59 x 10- 43 %. 

For yp — 1 we obtain in the whole interval of A, that is, yp G [0, oo), 

|pc(c) ^ (pp iff = 1; w = 50 °) 1 

I PD(C) A / e (A, A = 1; \J\P = loo) I 

Figj9] shows the percent difference (122]) for the entropy 
at | J\(3 = 500. For these values of yp and \J\(3, we obtain 

PD(S)il (A = 0.8, A = 0.51; \J\D = 50oi 

\PD(S) A / e (A > l,5i, 2 = 0.51; \J\0 = 50oi | 

The curve PD(S)p f x yp, at | J\/3 = 100 is plotted in Figj9jo in the interval yp G [0, 0.6]. 
The value of the entropy per site for both models, for some values of yp in the same figure, 
that is 

S§? G (~. 1,|7| =°,| 7 | =0-51, | 7 |=-2; 1^1/3 = 100^ ~ 1.037 xlO" 42 , (26a) 
=°>]7| =°- 51 »|7| =°; | ^= 100 ) ~ 1.203 xio- 42 , (26b) 


< 10~ 216 %, (24a) 

< 10' 42 %. (24b) 

per site, L = S, with A = 0.51 


-4.68 x 10- 70 %, 

(25a) 

io- 439 %. 

(25b) 
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and 


Sff G I -1, = 0.488, -- = 0.51, f- = -2; \J\j3 = 100 


D 


K 


\J\ 


\J\ 


\J\ 


S‘;r ( 1 , A = 0.488, j~r = 0.51,= 0 ;|J|/?= 100 


6.417 x 10” 2 , (26c) 


6.615 x 10” 2 . (26d) 


We also obtain that \PD{S)p^ ^yyy ~ 0.7, yyy = 0.51; \J\/3 = 500j < 10 8 and 

\PD(S)g > 1.51, ft = 0.51; | J\0 = lOo) < lO" 85 . 

Finally, for = 1, the percent difference (1221) of the entropy per site of the fer¬ 
romagnetic BEG model, with py = —2, and the spin-1 AF Ising model is such that 

I PD(S)g (i, ft = 1; W = 500) < 10-”*, and \PD(S)g (i, ft = 1; | J\/3 = loo) < 
10 -42 . Both inequalities are valid for p > 0. 


5 Conclusions 

We obtain the exact expression of the Helmholtz free energy (HFE) of the Blume-Emery- 
Griffiths (BEG) model in the presence of an external longitudinal magnetic field, for 
arbitrary values of the parameters in the Hamiltonian (J2J), valid for T > 0. The addition 
of the term — K(Sf) 2 (Sf +1 ) 2 to the Hamiltonian of the spin-1 Ising model with single-ion 
anisotropy term and external longitudinal magnetic field, enriches the phase diagram of 
the latter at T = 0[5j. Although our results are valid for the ferromagnetic (J < 0) and 
AF (J > 0) versions of the BEG model, in the present paper we restrict our discussion 
to the results of the ferromagnetic BEG model and take J = —1. From section 0 on, the 
parameters of the Hamiltonian (J2]) are scaled in units of \J\ and the inverse of temperature 
is measured as | J\/3. 

The thermodynamics of the ferromagnetic BEG model and of its respective phase 
diagram at T = 0, is analyzed for two different regimes of the y^-parameter: ^ < —1 
and > —1. In both cases, the value of the z-component of the magnetization, A4 Z , 
at low temperature, is determined by the ground states of the model presented in FigsJH 
The function -M-^) exhibits well-defined plateaux up to \ J\f3 ~ 1 for py > —1; however, 

for for py ~ —1 these plateaux are lost for | J\/3 ~ 10. 

The type of degeneracy of the ground state (exponential or non-exponential) along 
the lines that separate the phases in the phase diagram of the ferromagnetic BEG model, 
at T = 0, determines the behavior of the entropy and specific heat, both per site, of the 
ferromagnetic model at low temperature. We show that for y^y > — 1 the degeneracy of 
the ground state along the lines that separate the phases in Fig{T)o, at T = 0, does not 
diverge in the thermodynamic limit (N —> oo). In this region of values of yjy, the specific 
heat and the entropy, both per site, vanishes as T —y 0 in the presence of any external 
longitudinal magnetic field. 

The lines that separate the phases in the diagram of FigJTjr, where we have py < — 1 
and T = 0, are related to ground states with exponential degeneracy. We calculated the 
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degeneracy of the ground state along each of these lines and obtained their entropy per 
site at T = 0, and then compared these results with the ones derived from the HFE of 
the ferromagnetic BEG at very low temperature. We have a strong indication that the 
results agree in the limit of T —)■ 0. 

Finally for jjr = —2, the phase diagram of the ferromagnetic BEG model (see Fig), 
at T — 0, with ^ | is identical to the phase diagram of the spin-1 AF Ising model, 

with single-ion anisotropy term, at T = 0 (see FigjGjr). 

We compared the thermodynamics of both models at very low temperature (\J\/3 ~ 
500), when the largest contribution to their thermodynamic is expected to come from 
their respective ground states. At temperature \ J\(3 = 500, the function of both 

models coincide by at least 1 part in 10 14 . Such agreement decreases as the temperature 
increases. 

The ferromagnetic and AF nature of the models is not apparent when we measure 
Ai z in the interval of pj e [0, oo). They are distinguished, however, by the specific 
heat per site and the entropy per site for jyj ~ 0.51, when 2D ~ |J|, and for values of 
the external magnetic held such that 2 h ~ \J\, even at low temperatures, | Jf3 = 500. 
At this temperature and with this value of the parameter jjr, the specific heat per site 

under jjy = 0.51 are valued Cjff G ~ 10~ 213 and C 1 /™' 9 ~ 10~ 8 for the BEG and Ising 
models, respectively: their values differ by 205 orders of magnitude. We cannot explain 
this difference. When the value of the parameter yh- is such that 2D 3> | J|, the specific 
heat per site and the entropy per site are insensitive, at low temperature, if the model is 
ferromagnetic (J = —1) or AF (J = 1). 

The information about the exact thermodynamics of the staggered BEG model in the 
presence of an external longitudinal magnetic held can be obtained from the thermody¬ 
namics of the BEG model through the relation ([!)). The study of the AF BEG model is 
currently under way. We expect to present our results in the near future. 
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A The exact HFE of the D = 1 BEG model in the 
presence of a longitudinal magnetic field. 

The transfer matrix metliod|2. [3J 0] was used by Krinsky and Furman[TU] to calculate 
the HFE of the BEG model in the presence of a longitudinal magnetic held (L — 0 in 
Hamiltonian (1.1) of Ref. [10]). 

In this appendix we follow the steps of Ref. [5], but here we write the root of cubic 
equation (1A.4[1 with the largest modulus as a real number, obtaining the exact HFE of 
the model (]2]). 
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The partition function Z(J, h,D; (3) of the BEG Hamiltonian (J2]) is equal tojlj 


Z(J,h,D,K;f3) = Tr[V N ], 


(A.l) 


where N is the number of sites in the periodic chain and the matrix U for the symmetric 
Hamiltonian ([2]) is 


U (J,h,D,K;p) 


e -/3(J+2h+2D-K) e -/3(h+D) e ~/3(-J+2D-K) ' 

(3(h-\-D) ^ P{—h-\-D) 

e -p(-J+2D-K) e -p(-h+D) e -(3{J-2h+2D-K) 


(A.2) 


and [3 = in which k is the Boltzmann’s constant and T is the absolute temperature 
in kelvin. 

The matrix U (J,h, D, K] f3) is hermitian for any value of J, h, D, K and (3. Its 
three eigenvalues are real. The matrices U [see ea. (IA.2j) ] and T [in Ref.[10]] differ by 
a rearrangement of lines and the sign of the external magnetic field h. The partition 
function (1A. 1[) is an even function of h. 

In the thermodynamic limit (N —> oo), the partition function (I A. ID of the model and 
its HFE W are, respectively, 


Z(J, h, D, K] (3) = (Xi) N and W( J, h, D, K; /3) = —- ln[Ai(J, h, D; /?)], (A.3) 


for non-degenerate eigenvalues of U. We assume that Ai is the eigenvalue of matrix U 
with the largest modulus, root of the cubic equation 

- A 3 + PA 2 + QX + R = 0, (A.4) 

in which 

P = l + 2e~^ J+2D ~ K) cosh(2/3h) = fr[U] (A.5a) 

Q = Ae~ 2fiD cosh(2 f3h) sinh ~ + 2e~ m2D ~ K) sinh(2 f3J), 

(A.5b) 

R = -2e" 2/3(2D ' if) sinh(2 f3J) + ^-P^d-k) s in h(^j). (A.5c) 


The uniqueness of this eigenvalue for the matrix (lA.2|h is ensured by the Perron-Frobenius 
Theorem[nj for any temperature T. 

The roots of the cubic equation (jA.4[) are well known [IT?j. The root with the largest 
modulus is 


where 




cos(6 l ) = — -j^=^= 

V (~Q ) 3 


(A.6) 


(A.7a) 


18 





















with 


Q = 


3 Q + P 2 
9 


and 


R = 


9QP + 27R + 2P 3 
54 


(A-7b) 


The expression of Ai has cubic roots. Plotting the thermodynamic functions of the 
BEG model required numerical evaluation of that expression. In this work the CAS 
Maple has been used, and some spurious complex values (probably due to the way the 
cubic roots are handled by the system) appeared in the floating point evaluation even 
for 700 significant digits, so some caution had to be taken. The expression (1A.3|) of the 
HFE, valid at any temperature T and obtained from eqs. flA.6p ~ flA.7bp . is exact for the 
ferromagnetic (J < 0) and the AF (J > 0) BEG models in the presence of a longitudinal 
external magnetic field. 


B The states and energies of the dimers 

For spin-1 there are nine possible dimers in neighbouring sites of the chain. Those states 
and their corresponding energies, obtainable from eq.(j6]) are 


|D (A) ) i|i+1 = |0)j <8) |0)i+i and 
|-D (s) ) m+ i = |l)i <8> |l)i+i and 
\D {c] ) i)i+l = | - l)i <8> | - l)i+i and 


M) 

'i,i+1 

w 

SB) 
'M+l 

\J\ 

SC) 
" 2 , 2+1 

w 


= 0 , 


|-D (Ei) )m+i = |0)i <8> |l)i+i and 
\V (Ea) ki+i = |l)j O |0) i+ i and 

\D {Fl) )i, i+ i = |0)j ® | - l) i+ i and 
\D (F2) ) iyi+l = | - 1)* ® |0) i+ i and 

| J D (Gl) ) M+ i = |l)i <8) | - l) i+ i and 
|A> (G2) ) m+ i = | - l)i ® |l) i+ i and 


SB i) 
't,i+1 


Se 2 ) 


SBi) 
'i,i +1 


Sf 2 ) 
'i,i+1 


r (Gi) 
'i,i +1 

\J\ 

Sc 2 ) 

"2,2+1 

w 


2D 2 h 

K 

(B.l) 

+ PT“I7“ 

w 

(B-2) 

2D 2 h 

K 

(B.3) 

+ PT + PT 

Vy 

D h 

\j\ i jy 


(B.4a) 

D h 

i^i i jy 


(B.4b) 

n h 

= pi + pi’ 


(B.5a) 

n h 

z \j\ + w 


(B.5b) 

2D K 

■ + it ijt 


(B.6a) 

2D K 

■ + it Pi' 


(B.6b) 
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C The ground states and energies of the BEG model 
in the presence of a longitudinal magnetic field 

We assume that the chain has a even number of sites TV by letting TV = 2 M, in which M 
is a positive integer. In the thermodynamic limit (TV —> oo), we also have M —> oo. The 
ground state vectors at each phase in the phase diagram in Fig{Tt> at T = 0 are 


|^o)a = |0)i<8> |0> 2 <g> ■ ■ ■ <g> |0) 2M , (C.la) 

|'J / o)s = 11)i ® |1)2 ® • ■ ■ <8> |1)2M, (C.lb) 

|d / o)c = | — l)i | — 1)2 <S> • • • | — 1)2M, (C.lc) 

I x 3/o)^; = |0)i ® |1) 2 ® |0) 3 ® |1) 4 ® ® |0) 2 m-i ® |1) 2 m, (C.ld) 

^o)b' = |l)i <8> |0)2 <S> |1)3 ® |0)4 <g) • • • <g |1) 2 m-i ® |0) 2 m, (C.le) 


| v I / o)f — |0)i <£> | — 1)2 ® (0)3 <E) | — 1)4 <g) • • • <E> |0)2 m -1 <£> | — 1)2 m, (C.lf) 

^o)f' = I — l)l ® |0) 2 <8) I — 1)3 0 |0) 4 • • • <S> I — 1)2M—1 ® 10)2Af- (C.lg) 


One is reminded that S?\s)i = s|s)j, s G {0, ±1} and i G {1, 2, • • • 2 M}. 

The states |T 0 ) E >, \ and |^ 0 )g' have the same energies as the states |T 0 ) E , |T 0 )f 
and | v I / o)g, respectively. 

The values of the ground state energy of the phases A, B, E(E') and F(F' ) for ^ > 0, 
are, respectively, 


zdL 

\J\ 

El 

\J\ 

El 

\J\ 

H 

\J I 


0, 


2 M 

2 M 

2 M 


/ J 2D 2 h K \ 

VM + V\~V\~ \j\) ’ 

(D_ _ ±\ = E? 

VM \J\) \j\’ 


D h\ _ El 

\T\ + V\) ~ FT 


(C.2a) 

(C.2b) 

(C.2c) 

(C.2d) 


The phase diagram of the AF spin-1 Ising model, with a single-ion anisotropy term, 
in the presence of an external longitudinal magnetic field has an extra phase [5] at T = 0, 
the Neel state, given by the vector states 


|*0>G — 11) 1 ® | — 1)2 ® |1)3 ® | — 1)4 ® ® 11)2M—11 — 1)2M, 

d / o)G' = | — l)l ® |1)2 ® | — 1)3 ® |1)4 ® ® | — l)2M-l|l)2M- 


The ground state energy, in units of | J\, of theses states G and G 1 is 


\A 


Hi 

\J\ 


2M[ ~V\ + V\ 



El 
\J\ ' 


(C.3a) 

(C.3b) 


(C.4) 
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D Degeneracy of ground states along transition lines 
and on critical points of phase diagrams of the 
ferromagnetic BEG model with rfi < — 1 at T = 0 


Let us take the transition A ^ E/E' as an example for the calculation of the degeneracy 
of the ground states along the boundary of two phases in the phase diagram in FigOJ of 
the ferromagnetic BEG model at T — 0. The critical point T, excluded from this analysis, 
will be treated subsequently. For all points in the vicinity of (but not upon) this line, the 
chain may be found in either one of the following states, among those listed in Appendix 
□ 


\Vq)a = 

|0>! (8) |0) 2 ® |0) 2M , 


(D.l) 

|^o)b = 

|0)i<8) | 1) 2 <S» |0) 3 <8) | 1) 4 <S» ■ ■ 

■ • ® 10)2M—1 1 1 ) 2Mt 

(D.2) 

v Lo)e' = 

|l)i <8) |0) 2 <8> | 1) 3 <8> |0) 4 <8> ■ ■ 

■ ' ® 11) 2M—1 <8> 10 ) 2Af • 

(D.3) 


Hence, on the phase transition line, if the i-th site (for % G {1, 2,... 2 M — 1}, in which M 
is an integer) happens to be in the state |0);, the state of the (i + l)-th site may be either 
|0),: + i or |l)j+i. On the other hand, if the i-th site state is |l)j, the state of the (i + l)-th 
site must be |0)j+i. In other words, the energy constraint on this phase transition line 
imposes some sequencing rules on the states of the chain: each |0)j + i succeeds either a 
|0)j or a |l)j; and each |l) i+1 succeeds a |0)j. (Equivalently: |l) i+ i never succeeds |l)j.) 

Let and be the number of possible occurences of the one-site states |0) and 
|1) at the f-th site, respectively. The number of possible occurences for each state at the 
(i + l)-th site can thus be written as a recurrence relation of the form 

9 (i+ 1) = Tg®, (D.4) 


in which 


9 


(*) 


r w i 

3\o) 

, g w = 

' l ' 





, and T = 

± _L 

1 0 

»|i) J 


l 




(D.5) 


The configuration r/ 1 ) corresponds to the chain with one site (i = 1), in which we may 
have either the state |0)i or the state 11)i. 

Such recurrence generalizes to 


g(i+p) — JVg(i). 


(D.6) 


and hence we may write, relating the 1-st and (2M)-th (last) sites, 


g (2M) = T 2M-1(1) 


(D.7) 


The total number of states with the same ground state energy corresponds to the sum of 
all possibilities for the |0) and |1) states at the (2M)-th state. It is then equivalent to the 
L 1 -norm of the vector g^ 2M \ 


Qa^e/e' — \g {2M \ — S'lo) 


(2 M) 


+ g |i> , 


(D.8) 
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hence, we should turn our attention to the evaluation of the matrix power T 2AI 1 in (ID. 71) . 
The matrix T can be easily diagonalized, yielding 


A 


A^±E/E' 


;(l + \/5) 


0 


0 


-(1-^5) 


(D.9) 


and the corresponding matrix of eigenvectors, 

2 


D 


A^E/E' 


1 + ^5 1- V5 

1 1 


(D.lOa) 


so that 


We rewrite (ID. 71) as 


D A l±E/E' T D A^E/E' — h-A^E/E'- 


g{2M) = D A l ±E/E , K n a J e/e , D A^E/E'l^- 


(D.lOb) 


(D.ll) 


After some algebra, (ID. 8(1 yields the degeneracy of the ground states along the transition 
line A ^ E/E' in FigOi, at T = 0, 


n 


A^E/E' 


22 M 


5+3 ^h+v5r + w^h-v5' 2M 


10 


Thus the corresponding entropy for this transition is 

In Qa^e/e 1 


Sa^e/e 1 = bin 
M—>oo 


2 M 


= In 


10 


i + VE 


(D.12) 


= In A, 


0.48121, 


(D.13a) 

(D.13b) 


in which tp is the celebrated golden ratio. 

The reader should notice that the periodic spatial boundary condition on the chain 
(cf. section [1]) has not been used at all in the calculation of the degeneracy (ID. 12(1 and 
hence on the determination of the entropy per site (ID. 13bl) . We shall describe in what 
follows how the degeneracy can be calculated taking that condition into account; the value 
of the entropy per site (lD.13bjl . however, will not change. By identifying the (N + l)-th 
site of the chain with its 1-st site, the sequencing rules described in the beginning of this 
appendix should also apply to the IV-th and 1-st site. Here we have N = 2 M. Namely, 
if the chain has |l)i in its 1-st site, it cannot have |1 ) 2 m in its 2M-th site. On the other 
hand, if the chain begins with |0)i, there are no restrictions: it can end in either |0 ) 2 m c> r 
11) 2 m- The two branches of possibilities — chain starting with 10) 1 or chain starting with 
| l)i — can be described separately by the 1-st site conditions 


9 (1, '° = 


" 1 ' 


O 


. 9 (1W = 


1 

0 

_ 1 


1 


(D.14) 
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The version of (ID.71) upon boundary conditions reads 


= R 0 T 2M_1 gW’° + R 1 T 2M_1 g^\ 


(D.15) 


in which {i? 0 ) Ri} are matrices that implement the pertinent restrictions to each 1-st site 
condition. In the present situation, these matrices have the form 


Rn — 


1 0 
0 1 


and R\ = 


1 0 
0 0 


(D.16) 


The first term in (1D.15|) relates to the chain states with a |0)i state; Rq is simply the 
identity matrix, and there are no restrictions to the which state the 2M-th site may have. 
On the other hand, the second term in (ID. 151) relates to the chain states with a |l)i 
state; the effect of R\ on T 2AI ~ l g^’ 1 is that of discarding all possible chain states with 
a 11)2 m state, which would violate the energy condition/sequencing rules. (Notice that, 
without boundary conditions, we would have Rq and R\ both equal to the identity matrix, 
and (ID. 151) would be reduced to (ID. 71) . The 1-st site condition Eq. dD.5l) is equivalent to 
letting g (1) = g^' 0 + g^’ 1 -) The total number of states with the same ground state energy 
corresponds to 


O B.C. 

1 l A^E/E' 


' = 1 9 


(2 M) | 
B.C. Il 


(2M) 
BB.C., |l 


+ o (2M) 

^ Ub.c., |i); 


(D.17) 


analogous to dD.8() : here, g^c \o) an< ^ 9bc\i) are the components of g^c ■ We obtain 


(2 M) 


(2 M) 


o B.C. 
* l A^±E/E' 


22 M 


1 + 


2M 

5) + (1 


2 M 


(D.18) 


which differs from the degeneracy (ID. 121) : however, the same entropy per site is obtained 
from dD.181) as it is obtained from (ID. 1 21) . 


qB.C. 

d A^±E/E' 


lim 

M—>oo 


in n 


B.C. 

A^E/E' 


2 M 


Sa^±e/e', 


(D.19) 


given by dD.ldaf . 


The analysis for the critical point T is analogous. In the vicinity of (but not on) T, 
the available states are 


1^0 )a = 

= |0>i 

(8) 

|0) 2 <8) 

• • • ® |0) 2 m> 




(D.20a) 

^o)e = 

= |0)i 

® 

|1} 2 ® 

|0) 3 ® |i> 4 ® 

■ ® 

0)2 M- 

-1 ® |1)2 Mi 

(D.20b) 

'f'o)®' = 

= ll)l 

® 

|0> 2 ® 

|1> 3 ® |0} 4 ® ■■■ 

■ ® 

1/2M- 

-1 ® |0) 2 Af, 

(D.20c) 

h^o)F = 

= |0)i 

(8) 

1 _ 1)2 

® |0) 3 ® | - 1)4 

l ® 



(D.20d) 

T 0 )f' = 

= 1 - 

l)i 

<8> 0)2 

® | - 1 ) 3 ® |0)z 

1 <8> 

■■•® 

~ 1)2M-1 ® 0)2M- 

(D.20e) 


Hence, for the state T, the 7-tli site may be in any one of the states {|—1)*, |0);, |1);}, 
i & { 1,2,..., 2 M}. Moreover, for ie {1,2,..., 2 M — 1}, the state |—l)j+i can only be 
preceeded by |0)j, the state |0)j + i can be preceeded by |—1)* or |0); or |l)j, and the 
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state |l) i+1 can only be preceeded by |0)j in order to guarantee that the dimer has the 
least possible value of energy. (Equivalently: | —l) i+1 never succeeds either |—1)* or |l)jj 
|l),:_i_i never succeeds either |—l)j or |l) i .) The recurrence relation among the number of 
possibilities for each state at a site can be expressed as a recurrence relation of the same 
form as (|D.7|h but now with 



r (i) 1 

9\-i) 


i 

_i 

9® = 

(i) 

, i? (1) = 

i 




i 


The diagonalization of T leads us to 


and T = 


A r — 


2 0 0 

0-10 
0 0 0 


and 



0 1 0 
111 
0 1 0 


(D.21) 


(D.22a) 


(D.22b) 


By calculating g ^ 2A/ - ) 
at the critical point 


as in the previous situation, [see eq. f)D.12jl ]. we obtain the degeneracy 

T in FigHk, 


sij = l (2 2M+2 - 1), 

and the corresponding entropy 


Sj = 


lim 

M —^oo 


In ilj 
2 M 


= ln(2), 
w 0.69315. 


(D.23) 


(D.24a) 

(D.24b) 


Once again, the calculation of the degeneracy under the boundary conditions may 
be carried out in the same fashion as for the A ^ E/E' transition. Even though the 
degeneracy itself differs from that of the unconditioned case, 

= 2 2M + 1, (D.25) 


the entropy per site obtained is the same, 


S B.a 


In nf c - 

lim --— 

m- kx> 2 M 


Sj. 


(D.26) 


The degeneracies and entropies for the other phase transition lines in the phase dia¬ 
gram of the ferromagnetic BEG model with y/ < —1, at T = 0, and the critical point 1Z 
can be obtained in the same fashion. 
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It is important to point out that in order to calculate the entropy per site of the 
ferromagnetic BEG model along phase boundaries and multicritical points at T — 0 we 
can calculate the total number of degenerate ground states in the chain without taking 
into account the periodic spatial boundary condition. 
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Figure 1: The phase diagrams, at T = 0, of the ferromagnetic BEG model in the presence 
of an external magnetic field. In (a) we have ^ < — 1. The vector states corresponding 
to phases A, B, ■ ■ ■ , F/F' are given in appendix O The multicritical point 1Z is at 
jjt = 1 + In the phase diagram ( b ) we have yj| > —1. The phases A, B and C in this 
diagram are the same the ones that appear in diagram (a). The tricritical point in this 
diagram is at jjr — \ + ^ ■ 




\J\ \J\ 

Figure 2: The z-component of the magnetization A4 Z per site as a function of the external 
magnetic held |jy for the ferromagnetic BEG model. In (a), A4 Z is shown for = —2 
and several values of \ J\/3: 100 (solid line), 5 (dotted line) and 1.5 (dashed line). In (6) 
we have yjj = 0.25 with the curves plotted for distinct values of \J\/3: 1000 (solid line), 
100 (dotted line) and 10 (dashed line). In both figures we have J = — 1 and ^ = —1.5. 
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\J\ \J\ 

Figure 3: The z-component of the magnetization Ai z per site as a function of the external 
magnetic field py for the ferromagnetic BEG model in the region ^ > —1. In (a) we take 

|j| = —0.5. The curves are plotted for different values of | J\f3: 100 (solid line), 5 (dotted 
line) and 1.5 (dashed line). In ( b ) we have py = 1.8 and M. z is plotted for three values 
of \J\/3: 1000 (solid line), 100 (dotted line) and 10 (dashed line). In (a) and (6) we take 
J = — 1 and |jj = 1. 
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Figure 4: The entropy per site S as a function of the single-ion parameter is shown at 
| ,71/9 = 10 3 , with jj| = —1.5 and A = 0. We vary the parameter along the horizontal 
line py = 0 in phase diagram of Fig{l]a, in order to include in the curve the critical points, 
at T — 0, 1Z (jj| = —0.5) and T (py = 0). 
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Figure 5: The specific heat per site C of the ferromagnetic BEG model as a function of 
the external magnetic field |jy, with: J = —1, ^ = —1.5 and yj| = 0.25. In (a) and ( b ) we 
take |J|/3 = 10 3 and consider the variable in two intervals: [0.23,0.27] and [0.73,0.77], 
respectively. Off these intervals the specific heat per site is null at this temperature. In 
(c), the function C(A) is plotted at \J\(3 = 10. 




Figure 6: In (a), the phase diagram of the spin-1 AF Ising model, with single-ion 
anisotropy term, in the presence of a longitudinal external magnetic field[5], at T = 0. In 
(6), the phase diagram [TJi of the ferromagnetic BEG model, with = —2, in the presence 
of a longitudinal magnetic field at T = 0. 
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Figure 7: The difference (fl9|i between the ^-component of the magnetization, Ai z , as a 
function of the external magnetic field yyy, of the ferromagnetic BEG model with ^ = —2, 
and the spin-1 AF Ising model, with single-ion anisotropy term. In (a) that difference is 
shown for \J\/3 = 500; in (6), for \J\j3 = 100. 



\J\ 


Figure 8: The specific heat functions C per site of the ferromagnetic BEG model with 
yjy = —2 (solid line), and the spin-1 AF Ising model, with single-ion anisotropy term 

(dashed line), at \J\/3 = 100. Both curves are plotted for A = 0.51. 
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Figure 9: The percent difference of the entropy per site between the ferromagnetic BEG 
model with = —2, and the spin-1 AF Ising model, with single-ion anisotropy term, 

with t jt = 0.51. In (a), the curve is plotted at \J\j3 = 500; in (5), at | J\/3 = 100. 
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